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Abstract 



> 

\0 ■ We study the classical SU{N) (g> U{l)/SU{N - 1) U{1) Non Linear Sigma 

m : 

. model which is the continuous low energy effective field theory for N component 

' frustrated spin systems. The j3 functions for the two coupling constants of this 

m : 

0^ ■ model are calculated around two dimensions at two loop order in a low temperature 

' expansion. Our study is completed by a large analysis of the model. The (3 func- 

I . tions for the coupling constants and the mass gap are calculated in all dimensions 

^ ' between 2 and 4 at order As a main result we show that the standard pro- 

o ■ 

O ■ cedure at the basis of the expansion leads to results that partially contradict 

. ^ those of the weak coupling analysis. We finally present the procedure that reconciles 

^ ■ the weak coupling and large analysis, giving a consistent picture of the expected 
scaling of frustrated magnets. 
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1 INTRODUCTION 



The Renormalization Group (RG) approach to Critical Phenomena has provided one of 
the greatest success of Theoretical Physics in the past twenty yearsB . The most important 
concept that emerged is that of effective field theory: the recognition that, whatever 
complicated the microscopic physics of a system is, its low energy, long distance properties 
can be well described by a much simpler and universal theory. The concept of low energy 
effective theory is particularly transparent in Statistical Mechanics, through the celebrated 
idea of block spin procedure which allows to generate a long distance effective action 
by integrating out the short distance degrees of freedom of a microscopic spin system 
defined at the scale of some lattice size a = A~^. In practice, one has to have recourse 
to approximate methods among which the quantum field theoretical approach has been 
one of the most prolific and successful one'B' . In the latter, the block spin procedure is 
bypassed by a continuous limit in which is made the choice of the fluctuating fields together 
with operators satisfying the criterion of perturbative renormalizability. The resulting 
renormalizable action is then seen as the low energy effective theory, i.e. at scale p <^ A, 
of the microscopic theory since non renormalizable - or irrelevant - terms correspond to 
those that are discarded by taking the long distance limit A ^ ooS. Critical and, more 
generally, long distance quantities are then obtained from perturbative expansions of this 
low energy theory. As well known, the consistency of the perturbative approach relies 
on the existence of a fixed point, around which the strength of the interactions and, in 
particular, that of non renormalizable interactions can be carefully controlled. But the 
main task is then to ensure that the results obtained so far persist beyond perturbation 
theory. This generally consists in comparing results obtained from different perturbative 
approaches. A celebrated example for which such an approach has been particularly 
successful is the Heisenberg ferromagne 

However, many issues remain open in the context of more involved systems among 
which are frustrated spin systems (for a review see Indeed, the effect of competing 
interactions on the classical Heisenberg antiferromagnet has been the subject of a great 
interest, especially focused on the question of the universality class for its phase transition 
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in three dimensions. Previous studies based on the Non Linear Sigma (NLcx) model 
approach by Azaria et al. in as well as in the quantumliil' case 

have revealed the existence of a non trivial UV fixed point with 0(4) symmetry which 
governs the phase transition above two dimensions. But no stable fixed point has been 
found in the e = 4 — d expansion of the Landau-Ginzburg- Wilson effective action so that 
a first order phase transition is expected around four dimensions lli^ M . This situation is 
even more confused when one takes into account of experimental results that do not show 
any evidence of universality 0' S Q and those obtained with Monte Carlo simulations 
that are in favour of a continuous phase transition in d = 3 with exponents that differ 
from those of the 0(4) universality class0' IHI. Various scenarios have been proposed to 
account for these results but none of them is fully satisfactory. In any case, at the center 
of the debate is the existence of the 0(4) fixed point found in the NLo" model analysis. 
More precisely, the very question is to what extent can we believe on the results obtained 
with help of the NLcr model beyond perturbation theory. 

The validity of the low temperature expansion at the basis of the NLcx model approach 
has been intensively questioned. The possibility that massive or irrelevant operators as 
well as non perturbative configurations drastically affect the RG properties of frustrated 
and non frustrated systems is not excludedS' EH. In this work we would like to deal 
with the specific situation of frustrated spin systems. At the heart of the problem is the 
existence, in their RG phase diagram, of antagonistic fixed points that result from the 
competition between different kinds of order. The main consequence of the presence of 
several fixed points is that the status of some operator may change: being relevant with 
respect to a given fixed point, it may become irrelevant with respect to another one. In 
the low temperature expansion made around two dimensions, such an operator remains 
marginal so that one can describe this cross-over by means of a renormalizable theory 
interpolating between the different fixed points of the theory. However, when changing the 
perturbative scheme the power counting changes and there is no guarantee that irrelevant 
operators do not destabilize the 0(4) fixed point or, at least, reduce its basin of attraction 
so that the 0(4) scaling behavior becomes hardly observable. With this respect, one may 
question the existence of a universal theory that interpolates between the different fixed 
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points of the phase diagram. In order to investigate the rehabihty of the perturbative 
approach in such a multi-fixed point situation we shall, in the following, study extensively 
the NLcr model relevant for frustrated systems by means of low temperature and large N 

expansions. 

In the second section of this paper, we present the effective action that describes the low 
energy properties of frustrated spin systems which is the SU (N) (8) [/(!)/ SU {N — 1) <SiU (1) 
NL(7 model. This model differs from that used by Azaria et al. in their investigation of 
the low temperature properties of frustrated spin systems by the fact that the fundamen- 
tal fields, in terms of which the theory is written, span a representation of SU (N) rather 
than of 0{N). The reason we shall focus on this model will be explained in this section. 

In the third section, we proceed to an extensive analysis of the RG properties of the 
SU{N) (g) U{l)/SU{N - 1) (g) U{1) NLcr model in a double expansion in powers of the 
temperature T and e = d — 2. We give the P functions for the two coupling constants 
entering in this model, calculated at two loop order. The RG flow is solved and analyzed 
in 0? = 2 and d — 2 -\- e dimensions. We show that the (3 functions and the correlation 
length ^ interpolate between two fixed points, one with the CP^~^ symmetry, associated 
with high energy degrees of freedom, the other one with the 0{2N) symmetry, associated 
with the long distance behavior, which drives the finite temperature phase transition of 
frustrated spin systems in d > 2. As a main result, wc show that a new scale S, distinct 
from the correlation length, is generated in the RG flow. This scale, that characterizes 
the cross-over between low and high energy scaling behaviors, appears as the main feature 
of the low temperature RG analysis of frustrated models. 

It is the purpose of the fourth section of this work to address the question of the 
validity of the weak coupling analysis beyond the formal expansion in powers of T and e. 
This is why we study the SU{N) ® U{1)/SU{N - 1) U{1) NLcr model in an expansion 
in powers of 1/N . There appears here an important point, related to the way the cross- 
over length scale S takes place in this approach. Indeed, in the weak coupling analysis, 
this length scale is identified in the RG fiow analysis but is absent of the perturbative 
calculation which is completely free of any scale dependence. In contrast, in the 
expansion the situation is different and this length scale defines itself the way to perform 
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the perturbative expansion. As a result we show that, according to the value of S, one 
has to perform different expansions using different effective field theories in the vinicity 
of the two - CP^^^ and 0{2N) - fixed points. Then, the (5 functions for the coupling 
constants and the mass gap are calculated in all dimensions between 2 and 4 at order 
A careful comparison between the low temperature and large N expansions shows 
that there is an agreement of the RG quantities in the very neighborhood of the fixed 
points. In the last section of this paper we finally discuss the implications of our findings 
on the physics of frustrated spin systems. 



2 THE MASSIVE CP^-^ MODEL 

The most striking feature with frustrated systems is that long wavelength fluctuations 
are different, in nature, from the short distance ones. This is typically the case with the 
Heisenberg antiferromagnetic model on the triangular lattice. Starting with Heisenberg 
spins on the lattice, the strong correlations due to geometrical frustration induce a lo- 
cal order in which the spins adopt a 120" structure on each elementary cell. Whereas 
the fluctuating field at the scale of the lattice size a is a vector, it turns out to be, be- 
cause of frustration, a rotation matrix R of 5*0(3) at larger scales. It has been shown 



by Dombre and ReadlflJ that, as the consequence of the matrix character of the order 
parameter, the relevant action for frustrated Heisenberg models is a NLcr model defined 
on the homogeneous space 0(3) ® 0(2)/0(2). This model has been studied in the past 
with help of a low temperature expansion. As a main result of this study is the ex- 
istence of a stable fixed point in the neighborhood of two dimensions with increased 
symmetry 0(3) ® 0(3)/0(3) ~ 0(4)/0(3)i 0. Our aim being to investigate the re- 
liability of the low temperature expansion, we need a generalization of this model that 
allows a expansion. One of them is obtained when the action derived by Dom- 
bre and Read is written in terms of component real basic fields. This leads to the 
O(A^) (g)0(2)/0(A^-2) ®0(2) NLcx model which has been studied by Azaria et al.0 W 



Another option is to parametrize the long distance fluctuations with help of component 
complex fields of SU{N). Doing so we obtain the SU{N) ® U{1)/SU{N - 1) ® f/(l) NLa 
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model which we shall extensively study in this paper. The reason to favour this model 
is threefold: first, although we shall be interested in the classical case in the following, it 
appears that the SU{N) U{1)/SU{N — 1) ® U{1) NLcr model is more appropriate to 
investigate the quantum case. Indeed, there are reasons to believe that the low energy 
excitations in the disordered phase of quantum frustrated spin systems are spin one half 
deconfined excitations z called spinons01. On the other hand, this model has a more 
natural large expansion than the 0{N) CS) 0(2) /0{N — 2) eg) 0(2) version. Indeed, as it 
will be shown below, it exhibits the phenomenon of enlarge symmetry at the fixed point 
for any value of whereas this property is lost for the 0{N) ® 0{2)/0{N — 2) ® 0(2) 
NLcr model which has no special symmetry properties for > 3. Note finally that the 
SU{2) ® U{1)/U{1) and 0(3) ® 0(2)/0(2) NLa models, which are relevant for the case 
of Heisenberg spins, differ by their topological properties. Indeed, 5*0(3) having a non 
trivial first homotopy group, there exist point and line defects in two and three dimensions 
respectively which are not taken into account in the SU (2) formulation. It is not excluded 
that these topological defects drastically affect the physics of frustrated spin systems 



[25] 



We shall make however the hypothesis that it is not the case and concentrate ourselves 
on the multi-fixed point situation problem, leaving for future work the understanding of 
the role of topology. 

The partition function of the SU{N)^U{1)/ SU{N-1)®U{1) NLa model writes© © E^l 

Z = J Dz^Dz 6[zh - l) exp [-S] (1) 

with: 

/. [ \ d-2 \<i-2 / ^ \ 2" 

(2) 



S = j d'^x -^d^z^d^z + [z^ d^, zj 

Here z is a. N component complex vector submitted to the constraint z'^z = 1, A being the 
inverse of the lattice spacing a of the underlying spin system. The action (0) is invariant 
under the action of left SU{N) and right U{1) global transformations: 

z^U z e^^, U G SU{N) and 6 e IR . (3) 
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The ground state of @ is: 



(4) 



and is invariant under the action of a left SU{N — 1) subgroup and a diagonal U{1) group: 

(5) 



. / 2(iV-l) 



where -Dat-i is given by: 



I 1 ^ 



D 



N-l 



N{N -I) 



Action d^) thus describes the symmetry breaking pattern SU (N) ® [/(I) —>■ SU {N — 
U (1). It is a NLcr model defined on the homogeneous coset space SU (N) Cg) U{\)/ SU {N — 
1) ® U{1) that accounts for the physics of the resulting 2N — 1 interacting Goldstone 
modes. 

Let us describe the symmetries associated with action (Q) according to the values 
taken by the coupling constants /i and /2. When 4/i = /2, action (H) becomes gauge 
invariant with respect to the right U{1) local group and one recovers the CP^~^ model 
relevant for non frustrated magnets®. In all other cases the gauge invariance is broken 
as a consequence of frustration on the microscopic spin system. In the limit /2 ^ oo, the 
model becomes 0{2N) invariant as can be seen when parametrizing the z field in terms of 
a 2N component real vector field ^. The model thus interpolates smoothly, as one varies 
/2, between the two symmetric CP^'^ and 0{2N)/0{2N - 1) NLct models. 

The partition function (|l]) can be rewritten into another form which enlightens the 
symmetry properties and the role played by frustration in the model. As in the CP^~^ 



(6) 
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model, one can decouple the current-current interaction: 



UJ,, 



(7) 



with help of a Lagrange gauge field ^^0. After some standard manipulations, we obtain 
so a form that brings out more clearly gauge invariance properties: 



Z = J DzWz DA^ 6^z^z - exp [-S] 



with: 



where T>,, = d,, 



S 



iA^ and 



V,z^V,z + A'-'^A,' 



fi 



A' 

4/1 



(8) 



(9) 



(10) 



In the absence of the mass term M^, i.e. for 4/i = /2, action (|^) is gauge invariant 
with respect to the right U{1) local group and describes the pure CP^~^ model. When 
7^ 0, i.e. 4/1 7^ /2, action corresponds to a "massive CP^"^" model. We thus 
see that, in this parametrization, the effect of frustration is to give a mass to the gauge 
field A^. As an immediate consequence we can already anticipate that a new length scale 
has been generated by frustration. The next sections will be largely devoted to discuss 
in great details the physical and technical implications of the existence of this new length 
scale. 



3 WEAK COUPLING ANALYSIS 

3.1 Introduction 

The renormalizability of NLa models defined on coset space G/H around two dimensions 



has been studied by FriedanlH^). He has shown that the RG properties of a NLcr model only 
depend on the geometry of the coset G/H viewed as a metric manifold. In particular, the 
f3 functions for the coupling constants entering in the model are given, at two loop order, 
in terms of the Ricci and Riemann tensors on the manifold G/H so that renormalization 
depends only on the local properties of G/ H and is insensitive to its global structure. It is 



worth stressing that this weak couphng expansion is, in fact, a loop expansion ordered in 
powers of the temperature T. In symmetric models such as 0{N) or CP^~^ models, there 
is only one coupling which identifies with the temperature. However, in non symmetric 
models such as the massive CP^^^ model, there are at least two coupling constants. In 
this case, the resulting expansion is ordered in powers of an effective temperature but is 
not perturbative with respect to the ratio of the coupling constants. This is the reason 
why the RG functions may interpolate between fixed points that lie at finite distance. 
With this respect, the low temperature approach of the NLa model provides the natural 
framework to study frustrated systems whose phase diagram exhibits several fixed points. 
However, one has to keep in mind that the low temperature expansion misses terms of 
order e~^/"^ that include the global nature of the manifold - or of the massive modes of 
the theory - that may be important for the physics0. We shall return to this point in 
the next section when dealing with the expansion. 

Starting from the partition function (|1]) we can take advantage of the constraint of 
unit modulus to integrate out one degree of freedom of z. To this end, let us parametrize 
z as: 



(11) 



where tt is a A^ — 1 component complex vector field which, together with the scalar field vr, 
represent the 2N — 1 Goldstone modes. With this parametrization, the partition function 
can be rewritten in a more suggestive form: 



I Dn exp — f d^x Qijin) dn^diT^ 

J|7r|<l .27 



/|7r|<l I Z J 

which exhibits the geometrical nature of the theory. Here, vr* is a short notation for 
the 2N — 1 Goldstone modes. The action gijlir) dTc'^dir^ represents the line element on 
G/H equipped with the metric (7ij(vr). Thus, in the absence of the constraint |7r| < 1, 
action ([T2|) would describe a free theory on the manifold G/H. In the low temperature 
expansion this constraint is irrelevant at any finite order of perturbation theory since it 
gives contributions of order e~^^'^. Friedan has shown that the (3 functions which give the 
evolution of the metric with the scale A: 

A. = (13) 
9 



are given at two loop order by' 



1 „ 1 



= ^9ij - ^Rij - ^RT'Rjpqr (14) 

where Rij and Rjpqr are respectively the Ricci and Riemann tensors of the manifold G/ H 
equipped with the metric gij. As we shall now see, these intrinsic properties can be 
formulated in the language of group theory which enlightens the symmetry properties 
and provides a powerful framework for practical computations. 

3.2 Vielbein basis formulation 

The main advantage of the geometrical interpretation of the RG properties of NLo" models 
is that the Goldstone modes vr* (respectively the tangent vectors 5^vr*) can be viewed just 
as a particular coordinate frame on the manifold G/ H (respectively on the tangent space 
of G/H). For practical calculations, it is extremely convenient to work, in the tangent 
space, in the vielbein basis tu^ related to the d^ir^ via tu^ = cj/S^tt* such that the vielbein 
basis metric rju = gijlir) uj}uJj is free of any coordinate dependence. From the RG point 
of view, the main advantage of the vielbein formulation is that the geometrical quantities 
such as the Riemann tensor only depend on the Lie algebras Lie(G) and Lie{H), i.e. on 
the structure constants defined by the following commutation rules: 

[Ta, Tj] = fa/Tj , 

[Ta,n] = fab'T, , (15) 

[Tj,Tj] = fij''TK + fj/T,, 

where G Lie(if) and Tj G Lie{G)-Lie{H) are normalized according to Tr TjTj = —26ij. 
We shall now see how this works in the massive GP'^"^ model. 

Let us parametrize the coset SU{N)^U{1)/ SU{N — 1)®U{1) by choosing one unique 
element, i.e. fixing the gauge: 



(16) 



V 



_e-i7r e*^Vl - TrtTT 

This matrix is the same as for the coset space SU{N)/ SU{N — 1). This is the consequence 
of the fact that the coset space G ® X/ H ® X - where X is the maximal subgroup of G 
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that commutes with H - and G/H are topologically equivalent. In the following, we shall 
work directly in the coset SU{N) / SU{N — 1) keeping in mind that we look for an action 
which has SU{N) f/(l) as symmetry (isometry) group. 

The quantity L^^dL belongs to Lie(G') of G = SU{N) and we have: 



L-^dL = {L-'dL)G-H + {L-^dL)H 



(17) 



where {L ^dL)H is in Lie{H) oi H = SU{N — 1). One can now express the action (|]) in 
terms of the currents (ITTD as: 



S = -^ld''xTr[iL''dL)G-H] 



which is the general expression of NLcr models. It depends only on the currents in Lie{G)- 
Lie(if), the tangent space of G/H. This reflects the gauge invariance under the H group, 
i.e. the arbitrary in the choice of the ground state. Since L^^dL belongs to Lie(G), 
eq. ([T7|) can be rewritten as: 



L-'d,,L = LuiTj + niTa 



(19) 



where tu^ and are respectively the vielbein and the connection in the tangent space of 
G/H. They are given by: 

1, 



— Tr 
2 



— Tr 
2 



(20) 



Among all these currents, the U{1) one plays a crucial role. It is defined in terms of the 
z field as: 

1. ' 



X 



2N 



2"'~ 4ViV-l X^''^^'' 

while the currents u)^ in Lie(^f/(A^))-Lie(St/(A^ - 1))-Lie(f/(1)) verify: 

2" 



(21) 



(22) 



We can now write action (|]) in terms of these currents: 



S=-i (fx 



mE«)' + ^2«)' 



(23) 
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which defines the vielbein basis metric rjjj. Indeed we have: 



where 



with: 



1 

ViJ = -o 



S = \jd''x VIJ u^y^ , (24) 



r]^^I{T'T') + {r], - r7i)5,,5j^Tr(T-T^)] , (25) 
1 



(26) 



r/2 =2 



We recall that in our case T/ G Lie(S[/(A^))-Lie(5f/(iV - 1)), e Lie(t/(1)) and e 
Lie(5'f/(A^ — 1)) and that the corresponding algebra is given in (|1^). As already pointed 
out, whereas non frustrated models are defined with one coupling constant, the frustrated 
case is characterized by two independent coupling constants rji and 772- 

The crucial advantage of the vielbein basis formulation is that the Riemann tensor 
only depends on the structure constants /j/" of Lie(G) so that the (3 functions: 



= , (27) 



writes at two loop order: 



PiAv) = ^vij - ^Rij - A Ri'^'' RjpQR (28) 



where the Riemann tensor writes'0: 

1 



RlJKL — flj"'faKL + l^flJ UmKL + flMK — fKLKd) 



+ ^ iflKM + fhllK — IkMi) {^f/'b + fLJ^' — 

{fjKM + fhUK - fKMj) [fl^\ + JlI^'^ - f^\l) ■ (29) 

The indices a and {/, J . . .} refer to H and G ~ H respectively. The G — H indices are 
raised and lowered by means of rj^"^ and rjjj and repeated indices are summed over. 
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3.3 Results for the SU(N) U(1)/SU(N 



Ufl) model 



We are now in a position to study the RG properties of the massive CP^ ^ model. The 
computation of the Riemann tensor (eg . (|29|) ) leads to the following f3 functions: 



(30) 



with: 



V2 



Vi, 



N 



1 1 



1 iV 



^2 



^2 



AN -6 N —+ {3N — 



(31) 



??2 



where m = — rm- 

' 2{N - 1) ' 

In eq. (|30|) , the one and two loop order contributions to the {3 functions are ordered in 
powers of — = T, the effective temperature of the model. The long distance, low energy, 
physics is determined by the behavior of the RG flow in the A — > cxd limit while critical 
behavior is associated with the existence of UV stable fixed points of (POp. To exhibit the 
low T properties, it is convenient to make the following change of variables: 



T 



- = /i 



X 



4/1 

Vl /2 

In terms of these new coupling constants, eq.(^) reads: 

dT 



where: 



dx 



rp2 

(3t = —{N 

71 



1 +a;) + 



2^ 



(32) 



(33) 



- 1 + 2x + (3A^ - 1) 



(34) 



/3. 



T 



N - X {I- x) - N — X {I 



X 



IT 
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As usually, T is the relevant coupling constant that orders perturbation theory. The 
variable x measures the anisotropy of the model: as x varies one interpolates between 
the pure CP^^^ model, when x = 1, and the pure 0{2N) model, when x = 0. In the 
following, we shall study separately the d = 2 and d > 2 cases. 

3.3.1 Two dimensional case 

When e = 0, eq.(|33[) exhibits a line of UV fixed points at T = for any values of x. 
However they are not all true fixed points in the sense of limit of RG trajectories. To 
get the complete RG properties one has to solve the flow resulting from eg . (|33|) . In the 
following, we restrict ourselves to one loop accuracy since inclusion of two loops terms 
does not modify qualitatively our conclusions. 

The RG trajectories are parametrized by the RG invariant K^, given at one loop by: 

Kn = T, (35) 

1 — X 

in terms of which the solutions of eq.(^) are given by the implicit equation: 

K^lnX = An ix{l)) - An {x{\)) 

(36) 

= ^^^(A)i-i/iv 

with initial conditions x(l) and T(l) satisfying eq.(^). The function An{x) is given in 
terms of the hypergeometric function 2F1: 

^.W = ^-'-"'.f,(2.1-^;2-i;.) . (37) 

We show in Fig.l the corresponding flow diagram. 

In the UV regime, i.e. when A — > 0, all the flow lines with Kn 7^ reach the CP^^^ 
fixed point at x = 1 and T = which thus governs the short distance behavior of models 
with X 7^ 0. On the invariant line Kn = 0, corresponding to the 0{2N) model, the flow 
is governed in the UV by the 0{2N) fixed point at a; = and T = 0. In the infrared 
limit, when A increases, x{X) decreases to its asymptotic value at a; = 0. This means that 
all models with a;(l) 7^ 1 are asymptotically equivalent at long distance to the 0{2N) 
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model. We thus conclude that, while a priori we have a line of fixed points at T = 0, only 
two of them determine the RG properties. They are the CP^~^ and 0{2N) fixed points 
that govern the short and long distance behaviors respectively. The presence of these 
two competing fixed points is a characteristic feature of frustration where fiuctuations 
of different nature compete. We are thus led to expect that apart from the correlation 
length there should exist a new scale S, characterizing the cross-over between short 
and long wavelength associated with the CP^^^ and 0{2N) fixed points respectively. As 
a consequence, the scaling behavior of correlation functions will depend, in a non trivial 
way, on the relative magnitude of S and A^^. To see this, we have to go further in the 
fiow analysis. 

First of all, in order to perturbation theory remains valid one must have T(A) ^ 1. 
Given initial conditions on the fiow (x(l),T(l)) at scale A~^, this implies that A has to 
be much smaller than Aq defined by: 

An (x(1)) = Kn InAo . (38) 

This defines the correlation length = AqA"^ of the model: 

f - A-exp ( ^N{x{l)){l-xim 

^-^ e^p|^ x(i)i-v^r(i) j' ^^^^ 

which follows the standard RG equation: 

The general solution of eq.(^) writes in fact H = ^{Kj^) C, where ^{K]\f) is an arbitrary 
function of the fiow invariant Kf^f. This means that in a model with several coupling 
constants there may be others physical scales necessary to characterize the complete 
scaling behavior. We shall indeed see that if Aq defines the critical regime through the 
relation 1 ^ A ^ Aq it does not specify which fixed point governs the scaling behavior of 
the correlation functions. Let us now precise this point by defining three regimes according 
to the initial conditions on the fiow: 

• Given initial conditions very close to the CP^~^ fixed point with x(l) ~ 1, we define 
the region governed by the CP'^'^ fixed point to be such that x(A) is still close to 1. Since 
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TT 

in this reEfime Am(x) ~ — ; -, one finds that A has to be much smaller than 

^ ^ N{1 - x) 

We so obtain a new RG invariant scale S = AiA~^, solution of eq.(^), that defines the 
region governed by the CP^^^ fixed point: 1 ^ A ^ Ai or, in terms of the momentum 
p, <^ 'p <^ K. Notice that in this regime one has 1 — x(A) ~ T{\)/Kj^. 

• Independently of the initial conditions x(l), the 0{2N) regime is obtained for values 

of A ^ Ao such that x(A) ~ 0. One sees that it is sufficient to have A ^ Ai since An{x) ~ 

a;^^^^^ when x 0. The latter condition not only implies Km <^ 1 but also defines 

N-1 

the region governed by the 0{2N) fixed point: Ai ^ A ^ Aq or ^ p ^ 'Er^. 

• Finally, when Kn <t:l and x(l) ~ 1, one has ^ > S so that the 0{2N) and CP^~^ 
regimes can be observed when ^ p ^ and ^ p ^ A respectively. 

Let us resume our findings. We first found that the perturbation theory remains valid 
if one stays in a regime A <^ A,^ so that T(A) <^ 1. In symmetric NLcr models, such as 
when X = 1 or X = 0, the correlation functions display scaling behavior with logarithmic 
corrections that are governed by the unique non trivial fixed point of the theory. In non 
symmetric NLcx models, i.e. x(l) 7^ 0, 1, even when A ^ A,^, there remains a non trivial 
scale in the theory. In both regimes ^ J5 ^ and ^ p ^ A, i.e. when one 
lies near the 0{2N) or CP^^^ fixed points, the field theoretical approach to the RG is 
expected to give the correct asymptotic behavior of the correlations functions. 

The previous discussion can be formulated in terms of the gauge field formulation @) 
of the theory. This is most easily done in the limit 00, where S takes the nice form: 

S = (42) 

where is the gauge field mass (see eq.([TO|)). We see that ^ being fixed, the length S is 
associated with the energy scale where the effects of the U{1) current-current interaction 
become strong. We shall return to this point in more details when studying the 1/A^ ex- 
pansion. Let us just finally observe that the status of the f/(l) current-current interaction 
of eq.(|22D changes between the CP^~^ and 0{2N) fixed points. While it is marginally 
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relevant near the former, it becomes marginally irrelevant near the latter fixed point and 
is therefore a "dangerous irrelevant" variable. 

All the previous discussion is not qualitatively modified by inclusion of two loop order 
terms in the (3 functions. For completeness, we now give the expression of the correlation 
length ^ at two loop order: 

/ . _ ZN-A \ 

X 2iV2 



i — (,1-ioop G{x) 



3iV-4 



v 



(43) 



{l-x)^N{N-l)) 

where ^i-ioop is the one loop correlation length given by eq.(^), \E'o and \E'i being constants 
and G{x) being a pure function of x given by: 

'2F, (2, 1 - i; 2 - i; x) (1 - x)(3x + 1) N-1 + (2iV - 3): 



G{x) = exp / dx 



obtained by solving eq. 



2Nx 



2N{N -l)x{l~ x) 



(44) 



In the very neighborhood of 0{2N) fixed point, eq.(|43|) becomes to leading order in 



x: 



e = ^oi2N) (1 + ^) exp 



—X7l 



T{2N - 1){N - 1) ^ 



(45) 



T 



while in the neighborhood of the CP^ ^ fixed point, where 1 — x = — — h O(T^), one 



K 



finds to leading order in 1 — x and T: 

^ = ^cp^-i exp 



N 



(l-x)7r, T' 
-m — 



N^T 2tt 



(46) 



where Co{2N) and ^cpjv-i are the two loop correlation lengths of the 0{2N) and CP 
models0 §1 : 

( 1 

/ TT \ 



Af-1 



< 



^0{2N) — Ci^ 



^cpN-i = TN exp 



exp 



(A^ - 1)T 



(47) 



/ TT 



\NT 



where Ci^ and 6*2^ are constants that depend on the regularization scheme. 



17 



3.3.2 Above two dimensions 

In dimension d > 2, there is a pliase transition, as one varies T, from an ordered phase, 
with the breaking of the SU{N) ^U{1) symmetry, to a disordered phase. Indeed, eq. (|33|) 
admits, apart from the trivial T = fixed point, two non trivial fixed points with 

{T*pM-i = ^ - = 1) a'^d (T2(27v) = ~ 2(iV- 1)^ '^* " 

spectively CP^~^ and 0{2N) symmetric. In the infrared regime, the latter has one 
direction of instability and thus governs the long distance behavior while the former 
has two unstable directions and governs the short distance behavior. On the critical 
line connecting the two non trivial fixed points, the flow drives the system towards the 
0{2N) fixed point which is thus responsible of the phase transition of frustrated mod- 
els above two dimensions (see Fig. 2). The symmetry is increased at the stable IR fixed 
point for all values of A^. In particular when = 2, we find the 0(4)/0(3) universality 
class for frustrated Heisenberg spin systems in agreement with previous studies of the 



0{N) 0{2)/0{N - 2) 0(2) NLa model with = sHOl. 

As in two dimensions, there exists a specific scale = Ai^A"^ associated with the 
cross-over between the 0(2 A) and CP^^^ behavior. The scale can be again obtained 
by a careful analysis of the RG flow. As previously, we shall restrict ourselves to one loop 
accuracy. 

The one loop RG invariant flow takes now the form: 

l-l/AT 

Kn = - (T-r,(x)), (48) 

1 — X 



with: 



^^(^) = ]7rT(^-^)2^i(2'l-^;2-l;x) . (49) 



The critical line Tc(x) is obtained when Kn = 0, and Kn > label the lines in the 
disordered phase. With help of these quantities, the solutions of eq.(0) are given by: 

eAM{x{X)) + KN 



X- 



eAN{x{l)) + Kn 

(50) 
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where is given by eg . (^7\j . In the symmetric phase, when Kn > 0, the analysis of 
the RG flow follows straightforwardly from that oi d = 2. We obtain for the correlation 
length and the cross-over scale: 



K 



N 



e 1 + 



(51) 



NKm. 

which define, exactly as in = 2, the regimes of energy governed by either the CP^^^ 
or the 0{2N) fixed point. On the critical line, i.e. when Kn ^ 0, diverges but H 

has a finite limit: Sc = A^^ | ^ . This means that even on the critical line, 

the 0{2N) critical behavior of frustrated systems can only be seen asymptotically in the 
regime p <^ H^^. This fact may be of importance when studying frustrated spin systems 
on the lattice where one has to work with lattice sizes L ^ Sc to be able to observe 
0{2N) universal scaling. In the limit ^ cxd, Sc can be expressed in terms of of 
eq.® as: 

Let us finally give the expressions of the critical exponents: 

• At the 0{2N) fixed point, x is an irrelevant variable and scales to zero as: 



A 



-4'0{2N) 



x{\) ~ y—j , A > Ai, . (53) 

This defines the cross-over exponent 4>o(2N) which is obtained, at two loop, from eq.(^^: 

*o,.., = ^. + j^^.' + 0(.') . (54) 

The exponent i'o{2N) itself is, of course, given by that of the 0{2N) mo del0: 

^o(2JV) =e + ^^^ + 0(e3) . (55) 

• Near the CP^~^ fixed point, T and x are relevant. While the scaling exponent is 
given in the T direction by( 



UcpN-^ = e + ^e^ + 0{€') , (56) 
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we find that x{X) scales trivially: 

x{X) ~ 1 - (^-^^ , A « Ai, (57) 

so that (pcpN-i = e . This property holds at all orders of perturbation theory. This 
is easily seen on eq.(|3T|) where one can observe that the (3 function for f/2 orders with 
increasing powers of {fi2/fliY) p > 2, so that all its partial derivatives with respect to 
f]2 taken at 172 = vanishe. Hence the scaling of f/2 follows simply from dimensional 
analysis. This means that f/2 does not renormalize at sufficiently high energy and that 
the U{1) current has no anomalous dimension at the CP^'^ fixed point, a result which 
is completely consistent with the U{1) gauge invariance. 

3.4 Concluding remarks 

In this section, we have studied the massive CP^~^ model by means of a low temper- 
ature expansion around two dimensions. We have expressed the RG properties of this 
model in terms of local, geometrical, quantities associated with the manifold SU (N) 
U{1)/ SU{N — 1) ® U{1) viewed as a metric space. We have found that the f3 functions 
and the mass gap, calculated up to two loop order, interpolate smoothly between those 
of CP^~^ and 0{2N) models as one varies x. In addition, a careful analysis of the RG 
flow have revealed the existence of a new length scale S, characterizing the cross-over 
from short distance behavior, governed by the CP^~^ fixed point, to the long distance 
one governed by the 0{2N) fixed point. The appearance of this length results from the 
competition between fiuctuations of different nature and is reminiscent of the frustration 
of the lattice spin system. An important point is that, contrarily to the correlation length 
^, the length S is not related to a singular behavior of the RG flow so that perturbation 
theory is a priori well deflned in all the regime ^ p ^ A independently of the relative 
value of p and S, or equivalently, for any value of x between and 1. Indeed, the low 
temperature expansion which is at the basis of the NLcr model analysis is a loop expansion 
ordered in powers of T and e and is insensitive to the value of x. 

It is however worth to stress that all our results rely precisely on the perturbative 
renormalizability of the model in this double expansion in T and e. Under this hypothesis, 
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we expect that, at sufficiently low energy, i.e. when p -C A, the physics of the spin system 
is well described in terms of only two coupling constants x and T and is therefore universal. 
In this regime, the high energy degrees of freedom have decoupled from the low energy 
ones and all the remaining irrelevant operators give subdominant contributions to the 
scaling of physical quantities. The question that naturally arises is whether or not this 
universality property holds beyond perturbation theory and more precisely, to what extent 
can we control the effects of irrelevant operators beyond the double expansion in e and 
T. Indeed, one possible source of failure of perturbation theory in NLo" models is the 
neglected 0{e~^^^) terms that take into account of the global nature of the manifold. 
The best that one can do is to check the consistency of the low T expansion against other 
perturbatives approaches such as the expansion. It is the purpose of the next section 
to investigate the properties of the massive CP^^^ model through an expansion in the 
number of components of the fields and to study some aspects of the relationship between 
the large and weak coupling analysis of this model. 



4 THE LARGE N APPROACH. 



The 1/N expansion is a powerful tool on several accounts. On the one hand, it reveals 
by itself non perturbative informations about the model: mass generation, existence of a 
phase transition, etc. On the other hand, it provides a test of consistency of perturbative 
approaches, like weak coupling analysis. However, up to now, because of their strong 
analogies with gauge field theories in four dimensions and their simplicity, only symmetric 
models like CP^-i® B IH. IH' IH or 0(2iV)i M M © NLa models - with sometimes 
their supersymmetric counterpart - and Gross-Neveu mo del© m in two dimensions 
have been subject to a detailed analysis. In particular, it is just for the two last ones 
that the precise relationship between the large and weak coupling analysis has really 



been understood beyond leading order@. Such a relationship is still lacking in the case 
of more involved models like the massive CP^~^ model . 

As already emphasized, the main feature of this model is that the current-current 
operator {z"^ z)"^ which is relevant at the CP^^^ fixed point becomes irrelevant at the 
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0{2N) one. In the double expansion in T and e which is at the basis of the low temperature 
approach of the NLcr model around two dimensions, such a term stays marginal and is 
consequently renormalizable. This is why the low T expansion allows to interpolate 
between the two fixed points of the theory. However this situation is particular to the 
weak coupling expansion and one does not expect that this will be the case in other 
perturbative scheme. Indeed, we shall see that the standard expansion is unable to 
describe the whole phase diagram of the massive CP^~^ model. The reason for this is 
that the cross-over scale that emerged from the RG flow analysis of the weak coupling 
expansion enters in the expansion in such a way that it provides a boundary to the 
domain of validity of the standard large perturbative approach. We shall show however 
that by carefully taking into account of the presence of this new length scale, it is possible 
to capture the correct renormalization of the parameters of the massive CP^^^ model in 
the large analysis at least around each of the two flxed points of the model. 

Let us flrst explicit the dependence of the coupling constants /i and /2. The action 

writes now: 

AT „ .._o AT / " ^2l 



(5^ 



where the z is a A^ component complex fleld still subject to the constraint z'^z = 1. The 
starting point of the large A^ expansion is to take into account of the latter constraint 
with help of a Lagrange multiplier fleld A with a resulting effective action: 



S 



(fx 



d^z^df^z 



iXiz^z-A'^-^—] + — A 



N 



2t 



xt 



2N 



da Z 



(59) 



with t = fi and x = 4/1//2. In the absence of the current-current interaction, quartic in z, 
one could have integrated exactly over the z fleld to obtain an effective, non local, action 
for the A fleld from which results the well known 1/A^ expansion of the 0{2N) model. The 
difficulty with the massive CP^~^ model comes precisely from the U{1) current-current 
interaction. The standard prescription for dealing with such a term is to introduce a 
Lagrange fleld A^j, as in the pure CP mo In this case, the resulting action 

is given by: 



S = j (fx z'^ (-Vl - a) z + A 



2 + 2t 



(60) 
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where = {l — x)/xt and T>^ = d^j_ + iA^ is the covariant derivative introduced in 
section(2). Now the integration over the z field can be exactly done and the resulting 
partition function writes in terms of A and fields as: 

Z = j VXVAf, e~^^o (61) 

with: 

So = Trln(-Pj - ^A) + A'^"'- d^'x X + A^-^— / d^'x A^^ . (62) 



The set of eqs.(|6Ti52) have been obtained first by Campostrini and RossiEU who 
have studied the resulting expansion in dimension d = 2. They have computed 
the mass gap and the f3 functions for the coupling constants t and x at leading order in 
As a result, they argued that the mass term does not renormalize, a result 
which contradicts our two loop expressions of eg . (|34D . The origin of this contradiction 
between the weak coupling and large N analysis lies in the fact that the decoupling of 
the current-current interaction is not an innocent procedure. In particular, the use of the 
field theoretical approach on the resulting field theory ( |62D to obtain the correct evolution 
of the two coupling constants makes sense only in a regime of energy much greater than 
any other mass scale of the theory, i.e. in the vinicity of the CP^^^ fixed point. This 
fact was missed by Campostrini and Rossi. As a consequence, their results for the mass 
gap or the (3 functions are not correct in the large mass limit. In this case indeed, 
one has to return to the form (|S3) of the action and considerate the current-current term 



xt ( ^ \ ^ 

z\ as an irrelevant operator in the neighborhood of the 0{2N) fixed 

point. 



In the following, we shall study action ( |5BD in the large N expansion in = 2 and 
d > 2. Some of our results in ci = 2 have already been obtained by Campostrini and 
Rossill3. To be self contained, we shall present them together with new ones. 

4.1 The iV ^ oo limit 



As readily seen on eq.(^), the theory is exactly soluble in the limit N oo since the 



functional integral is dominated by the saddle point: 
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r 5So 

5X 
SSf) 



(63) 



. 



Searching for solutions with A^^ constant one finds A^ = so that the saddle point is 
the same as for 0{2N) and CP^-\ Eq.(|3D then leads to: 

d'^k 1 A'^-^ 



{27rY A;2 + m2 2t 



(64) 



with A = im"^ 



In dimension d > 2, the integral (|6^) is UV divergent, and we shall use a cut-off 
regularization procedure. One has to consider separately the d = 2 and d > 2 cases. 
• In d = 2, eq. (|6^ gives the mass gap equation: 

A 



TT 



In- 



t m 

from which follows the /3 function for t to leading order (A^ = oo): 



(65) 



vr 



(66) 



which expresses the property of asymptotic freedom in two dimensions* 
• In d > 2, eq.(p^ writes: 



d^ 

r 1 - - 

9 d-2 



(477)^/2 

2 



(67) 



^ - 2 

with tc = — r— and Kd = , . ^ ./n^/ , / n • 
2Kd (47r)'^/2r(rf/2) 

One can then deduce the (5 function for t: 

P, = -(d-2)t(l-^y (68) 

The /3 function for x follows trivially from the fact that does not renormalize at this 
level since A^ does not participate at the saddle point. One thus have: 



Pm^ = {d- 2)M' 



(69) 
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The set of eqs.( |6q - |69|) gives a complete description of the RG properties of the model 
at = oo and is completely consistent with those obtained in the weak coupling analysis 
of the previous section: 

• In dimension d = 2, the expressions of the mass gap and the (3 function for t agree 

with those obtained from the weak coupling expanded at leading order in A^. Moreover, the 

non renormalization of the mass term at this order is consistent with the fact that 

X — X 

identifies with the weak coupling one loop flow invariant at = oo since = = 

xt 

K^' (see eq.(il)). 

• In dimension d > 2, there exist two non trivial fixed points with t = and = 
or = oo. The fixed point with = oo governs the phase transition for all models 
with M V since is a relevant variable. The critical exponent u ^ = d — 2 is trivially 
the same as for the 0{2N) and CP^~^ models while the cross-over exponent = rf — 2 is 
just given by dimensional analysis. Both results agree with those of eqs. (|5^p7D at leading 
order in e and A^. 

As seen, there is a complete agreement between weak coupling and large A^ approach 
to leading order in 1/A^ and e. However, the N = oo limit is special since it does not 
involve the U{1) current-current interaction which is crucial in the CP^^^ model. This 
interaction will express itself as soon as 1/A^ corrections to the saddle point will be taken 
into account and will reveal the non trivial structure of the RG properties of the model. 

4.2 The 1/N corrections: general discussion 

Let us now investigate the effects of fluctuations around the saddle point. The Feynman 
rules for the 1/A^ expansion of action ( |60D follow from those of the pure CP^~^ model 
(see for instance Ref. ||31||). They are summerized in Fig.3. 

The inverse propagator of the A field is given as in the CP^^^ model by: 

^^^^^ " / {2^Y (g2 + rn?)(^(^q + kf + w?) ' ^^^^ 
The only difference with the CP^^^ model comes from the inverse gauge field propagator 
which acquires a longitudinal part: 

r..(fc) = [5,. - ^) T^{k) + A'-'M^'^ (71) 
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with: 



4 r(2-rf/2) 



m 



,d-2 



+ 



Am/ 



(72) 



d-l (47r)C?/2 ' d-l 

With help of the rules of Fig.3, the corrections can be derived in a standard way 
for all Green's functions of the theory with the restriction that no internal spinon loop is 
allowed since they are already taken into account in the definition of the A propagator. 
Before turning to the practical computations of the corrections to the mass gap m 
and to the (3 functions for t and or it is worth at this stage to precise from eg. ([72|) 
what is the expected domain of validity of the decoupling procedure. To see this, let us 
consider the short distance behavior of the inverse gauge field propagator involved in the 
computation of the correlation functions. In the large k limit, one has: 



with A, 



T^{k) ~ A' 

fc— >oo 



d-2 



m 



d-2 



and Sd 



4 r(2-rf/2) 
d-2 (47r)^/2 

r(rf-l)sin [7r(rf-2)/2] 



in d 



in d> 2 



(73) 



2(d - l)5rf ^" 2Ti{T{d/2)f 
It appears from expression ( [73D that the scaling behavior of the coupling constants 



(x, t) can be extracted from the UV behavior of the field theory ([6^ ) when: 

in d = 2 



1 

« -^In^ 

27r m"^ 



(74) 



< Ad in > 2 . 

This defines the energy regime (I) where the decoupling procedure is meaningful: 



p ^ m exp (ttM^ 



p > A 



M 



2\ 



m 



m 



d 



d>2 . 



(75) 



Expanding the second equation to leading order in e, one sees that the bounds defined 
in eq.([75D identify with those obtained in the weak coupling analysis - eqs. ( ^p2D - that 
define the regime governed by the CP^^^ fixed point. Using the saddle point equations 
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(see eqs. (|65|j67D ), one can translate eq. (|74[) in terms of t and as: 

tM^ < 1 in ci = 2 



d-2 



(76) 



in d> 2 



4(d-l)5d 

There exists a second regime (II) where it is possible to extract the renormalization 
of t and X. This is the regime governed by the 0{2N) fixed point defined by: 

p <^ m exp (vrM^) in d = 2 

(77) 



p < A 



2\ l/d-2 



in d>2 



which, in terms of t and M^, writes: 



tM^ > 1 



in d = 2 



d-2 



m 



> 2 . 



4(d-l)5d 

In this regime the field does not propagate anymore and the decoupling ( |60D is of no 
use to determine the renormalization of t and x. One has to return to the original form 
of the action (|59D. Near the 0{2N) fixed point, the current-current U{1) interaction is 
irrelevant and can be considered as a composite operator in the 0{2N) theory so that 
its effect on the physics can be computed perturbatively in a double expansion in x and 



In the following, we shall consider perturbation theory in the regimes (I) and (II) 
separately. We shall compute the P functions for the coupling constants x and t and the 
mass gap m at leading order in in = 2 and d > 2. Comparison with the results 
obtained from the weak coupling analysis of the preceeding section of this work will also 
be made. 



4.3 The CP 



N-l 



regime 



This is the regime (I) defined by eg. ([7^) . As already quoted above, perturbation theory, 
as given by the Feynman rules of Fig.3, is well defined. We shall assume in the following 
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that the theory is renormahzable so that only renormahzation of spinon and gauge field 
masses and fields are needed. We shall thus concentrate on the two point functions of the 
spinon and gauge field at order 

(79) 

[ G^Zi/NiP) = (A(O)A(P)) . 
The renormahzation procedure is standard. As usual, we shall consider the leading 
order of the Taylor expansion near p = of Gfj^{p) and which suffer from 

divergent corrections which will be eliminated by appropriate redefinitions of the masses 
and fields. 



4.3.1 The two point spinon function 

The diagrams that contribute to the correction of G^y^ip) are given in Fig.40. The 
self energy is defined by: 



i/N f/ p2 _|_ ^2 N + m? "^^'p^ + m? 
or in terms of the 1-P.I. function: 



^0) 



Tfl{p) = P' + + -nP) ■ 



N 



After lengthly but straightforward calculations we have obtained for S(p): 

S(p) = Ap^ + Bm^ + 0{p^) 



12) 



with: 



A 



A-d f d'^k Tl\k) 



4(d-l) f d'^k Tl\k) A 



{27rY {k^ + m^f 



2-d 



(27r)^A;2 + m2 J {27iY k^ + m 



d'^k 



B = 3 



d'^k T^\k) 
{27tY + m^f 



d^k 



1 



{2ttY A;2 + m2 



+ (d-l) [A - {d - 1)A''-'M 



,-2,,2^x\0)\ f d'^k Tl\k) 



{2tiY P + m2 
(83) 
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Defining: 



with the following prescriptions: 



rg)(/ = o) = m?, 



^4) 



dV%\p) 



p=0 



we obtain: 



A 
N 



1 + 



^6) 



B-A 

N 



We shall now consider the cases d = 2 and d > 2 separately. 
Two dimensional case. 
In d = 2, we find for Z and m|.: 



1 



1 

N 



A2 A2 

In In — - \n{2nM^ - 2 + In 



1 , A2 
In — - 



AtxM^ 



17) 



rrij^ = m 



1 / 

1 + — In In — + (3 - 27iM^) \n(27iM'^ - 2 + In — 1 



These expressions agree with those of Ref . [p7[] . Using the saddle point equation (|65D, one 
can express the mass gap m|j as a function of t and M^: 



1 + 1 (^In ^ + (3 - 27rM2) ln(27rM2 - 2 + 



(88) 



Using the N = oo result (|66D and the fact that Pj^p = 0{1/N) since does not 
renormalize at the saddle point one obtains the (3 function for t: 



TC 



1 + 



t 



2-kN 



1 



3 - 27rM2 
1 + t(M2 - l/7r)_ 



^9) 



This result has been first obtained by Campostrini and Rossi0 
• Above two dimensions. 

Let us begin by giving the explicit expressions of the field and mass renormalizations 
in = 3: 
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= m 



32 A 16 A 256 /l vrAM^X , /A , m 

■ In — + + —— In — + M^A - 



37r2A^ m ttN m -k^N V3 
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which leads to the (3 function for t: 



A = - 



4 



M2 



^ 7r2 1^^ ' A^l/16 + M2 + l/t-l/7r2^ 



1 + 



16 



3A^7r2 3A^7r2l/16 + M2 + l/t-l/7r2 
In d > 2, we can obtain the /? function for t from eqs. (|67D and (^31): 



2tx 



(90) 



(91) 



t - 



(P-d-2 {d-l 



M2 



iV 



Arf + M2 + l/t-l/t, 



_ 2 ^ 4(rf-l)(rf-2)g, ^ 2(^2 _ i)(rf - 2)ir, 



Nd 



Nd 



^d + M2 + l/t- l/t,_ 



(92) 



4.3.2 Gauge field renormalization 

We shall now investigate the role of fluctuations on the gauge field at order 1/A^ in two 
and three dimensions. This will provide the renormalization of the gauge field mass M"^ 
from which will follow the renormalization of the complete set of coupling constants of 
action (|5^. From the renormalization point of view, the situation is formally equivalent 
to what happens in Quantum Electrodynamics. In particular, Ward identities stated in 
the treatment of gauge field with a non- vanishing mass for the photon ensure the absence 
of a counter-term for the longitudinal part of the gauge field propagator in such a way 
that the renormalization constant for the gauge field A^^ Zjs,^ and for its mass, Zm'^ verify 
Z^a^Zm'^ = 1- It follows from this that the inverse gauge field propagator at order 1/A^ 
can be written: 



(93) 
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where r±i/N{p) is the transversal part of the inverse gauge field propagator at order 
We define the renormalized inverse gauge field propagator by: 



with: 



where Za is determined by the prescriptions: 



r±ij(p = 0) = M' 



(94) 



(95) 



R 



(96) 



and: 



dV 



p=0 



127rm|. 



1 



in = 2 



in d = 3 . 



(97) 



The diagrams contributing to the renormalization of the gauge field at this order 



can be found in Refs. p7|, |4^ , |45[| and they are reproduced in Fig. 5. The renormalization 
of T x^{p) is provided by a Taylor expansion in power of the external moment p of the 
diagrams. We then need the expressions of the transversal part of the bare inverse gauge 
field propagator at small p\ 

^2 



r±(p) 



r±(p) 



p- 



+ + 0(p^ 



p" 



+ + 0(p^ 



md = 2 



in (i = 3 



(9^ 



p^o 247rm 

which, taking into account of the spinon mass renormalization at order given by 
eg. (PB|) , leads to: 



r±(p) 



r±(p) 



+ .A o + + O(p^) 



P^ B - A , ^/ 4x 



md = 2 



in d = 3 



(99) 
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The calculation of the counter-terms for the inverse gauge field propagator at order 
is very cumbersome and details are given in the Appendix. 
We have obtained the following expressions: 



IZTTmp, N 



r±i/iv P = p - — r— m d 



(100) 



which, together with eg . (|99|) , show that the counter-terms associated with the transversal 
part of the inverse gauge field propagator are just cancelled by the renormalization of the 
spinon mass in such a way that the gauge field renormalization is trivial: 

Za, = 1. (101) 

We thus find that the gauge field mass does not renormalize at order 1/A^ in c? = 2 and 
d = 3. This result probably holds for any dimension between 2 and 4. However, it is worth 
stressing that this result is true only in the small tM'^ limit where the above calculation 
makes sense. The non renormalization of at order 1/A^ and the relation between 
and X - = {1 — x)/xt - induce a non trivial renormalization for the variable x given 
by: 

4.3.3 Comparison between large and weak coupling analysis 

At first sight, the /3 functions for t and x seem to interpolate between the CP^^^ and 
0{2N) models: they behave smoothly as functions of tM"^. The f3 function (|89|) for t 
allows to recover that of the CP^^^ and 0(2A^) models as goes to zero or infinity 
respectively. However, as discussed above, the whole calculation and, in particular, that 
involving gauge field loops is valid in a narrow range of energy defined by tM ^ ^ 1 which 
also implies 1 — x ^ 1. We shall now see that it is only in that case that the results of the 
1/A^ approach agree with those of the weak coupling analysis. In contrast, we shall show 
that the naive extrapolation of the results for large tM"^ disagrees with the weak coupling 
ones. 
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Let us first show that there is, in the CP^~^ regime, a good agreement with the weak 
couphng results of the preceeding section. We consider now separately the cases d = 2 
and d > 2. 

Let us consider the f3 functions (0) obtained in the weak coupling analysis in two 
dimensions. These equations can be expanded in powers of by doing the substitution 
T —>■ t/N. The resulting expressions can then be expanded in powers of 1 — a;. Keeping 
only terms of first order in this parameter one obtains the following expressions: 



t 2t^ 
Px = — (1 - a;) - — -(1 -x) 

77 77^ A/ 



(103) 



Now one can formally expand the f3 functions (|89D and (|102|) obtained in the 1/A^ approach 
in powers of t up to order t^. Then is replaced by its expression as a function of x 
and t. Keeping again just terms of order 1 — x, we recover the set of equations ([1031). 



This matching between the large and weak coupling expansions also works for the 
mass gap. The expression (^) of the mass gap obtained in the large N analysis can be 
expanded in powers of t and re-exponentiated in powers of 1/A^. This leads to: 

« = CA- Jexp(^)exp (-(i^l„±) (104) 

which coincides at leading order in t with the expression (|46|) of ^ obtained in the weak 
coupling analysis in the CP^^^ regime. 

In dimension greater than two, one can compare the exponents VcpN-i and (pcpN-i 
that govern the leading scaling behavior of x and t at the CP^~^ fixed point. The 
exponent v^pN-i obtained in the 1/A^ expansion from eq.(p2])@: 



coincides with that obtained in the weak coupling analysis (eq.(^)) at order 1/A^ and 
e^. The cross-over exponent (pcp^-i can be obtained in the large N limit from eq.( p.02| ) - 
(j)(jpN-i = d — 2 - m agreement with the weak coupling result. We can therefore conclude 
that, in the CP^~^ regime, the large N and weak coupling expansions coincide. 
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Such an agreement does not exist in the 0{2N) regime. In two dimensions, let us 
consider the /3 functions obtained in the weak couphng analysis and re-expanded in powers 
of 1/N: 



(106) 



xt 



1 + 



t 



TlN 



These expressions have to be compared to the f] functions obtained in the large N analysis 
expanded in powers of l/tM"^ since the condition tA'P ^ 1 defines the 0{2N) regime: 



TT 



1 



X 



A = - 1 - T7 + 



13. 



xt 

TT 



N N 



1 

N 



2n^N 



'1 



X 



(107) 



XV 



27r^N 



As readily seen, the two set of equations ( |106| ) and (|107|) differ for the /? function for x 
even at first order in t. This disagreement persists, of course, for the mass gap by a term 
which is linear in x/N. 

In dimension greater than two, let us compare the critical exponents obtained from 
the (3 functions of eqs.(p2|, |102| ) with those obtained in the weak coupling expansion. 

The critical exponent i'o{2N) as obtained from eq.(^ 



d 



4(rf-l)gd 
Nd 



(108) 



coincides with that obtained in the weak coupling expansion (eq.(^5D) at order 1/A^ and 
e^. However, since the gauge field mass does not renormalize at order 1/A^, the cross- 
over exponent is given by dimensional analysis: 4>o{2N) = d — 2 + 0{1/N'^), and thus 
disagrees with the weak coupling expression of eq. (^4]) . There is, in the 0{2N) regime, a 
disagreement between the 1/A^ expansion performed with help of the decoupling procedure 
and the weak coupling expansion. The reason is that the decoupling procedure does not 
take properly into account the dangerously irrelevant character of the current- current 
operator. The agreement we have obtained for the critical exponent i'o{2N) stems from 
the fact that the exponent u is determined by the asymptotic scaling in the relevant 
direction t at x = 0. 
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Naturally, this disagreement was expected since as previously discussed, the decoupling 
procedure is valid provided tM^ ^ 1. We shall now see that, provided one performs a 
suitable expansion in the vinicity of the 0{2N) fixed point - i.e. when tM'^ ^ 1-, 
the weak coupling and the large analysis are consistent. 



4.4 The 0(2N) regime 

In this regime, the gauge field no longer propagates and perturbation theory defined 
by the Feynman rules of Fig. 3 is meaningless. In particular, the decoupling of the current- 
current term (^z'^ dfi is of no use for exploring this low energy regime and one has to 
come back to the original form of action (|59|). Simple power counting around the 0{2N) 
theory indicates that (^z^ has scaling dimension 2 — d and is therefore an irrelevant 

(resp. marginal irrelevant) operator in dimension d > 2 (resp. in d = 2). Its effects on 
the low energy physics can however be computed in a double expansion in x and 
around the 0{2N) theory as soon as x ^ 1 and ^ lE3 EB. In the following, we shall 
compute the corrections to the RG functions at first order in x and The Feynman 

rules for the double expansion follow from those of the expansion of the 0{2N) sigma 
model. Propagators of the spinon z and A fields as well as the Xzz vertex are unchanged. 
The only modification is the presence of a new four point vertex F^j which has been first 
introduced by Chubukov et al.' — ' ED in their study of the quantum version of the model 
in d = 2 + 1 (see Fig.6): 

F^ ih, h, h, h) = A'-''— [h + h)^ {k2 + h)^ 4 (fci) 4 (fca) ih) zp ih) . (109) 

This operator being irrelevant we expect that higher order derivatives irrelevant operators 
will be generated by renormalization. However as we are only interested in the leading 
correction in x we can omit these possible additional terms. 

As in the preceeding subsection, we shall consider the renormalization of the spinon two 
point function and the four point vertex function F^^"^) at first order in x and 1/A^ from 
which the RG functions will be computed and finally compared with those previously 
obtained in the weak coupling analysis. 
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4.4.1 The two point spinon function 



(2) 

The diagrams contributing to the two point function G\L{p) at first order in x and 



are depicted in Fig. 7. The coefficients A and B entering in the expression of the self 
energy of eq.(]8^) are now given by: 



A 



4-d f d'^k Tl\k) 



d 



2^2 



= 3 



dn-. VI' (k) 



{3-d) 



d'^k 



{2nY F + m2 



^ (0) 



d'^k 



(110) 



In terms of this quantities, the field and mass renormalizations are given by eq.( 
• In two dimensions, we find: 



Z = 1 



N 



1 , , A2 , A2 
- In In — In — - 

2 An 



'1111 



niji = m 



1 A2 1 A2 1 xt fA^' 
1 + TT In In — r - — In — - + —— In — - 



N 



m? N m? N An 



With use of the saddle point equation (|65|) , the mass gap m\ writes in terms of x and t 
as: 

9 9 r 1 / 27r 27r X7r\l , , 

™?, = m^[l + -(ln--- + -)J (112) 

from which follows the j3 function for t: 



n\ n)^ 2n^N 



(113) 



• In dimension greater than two, the divergent part of A and B are easily computed 
and, omitting the details, we find for (3t: 



t 



Nd 



Ntr. 



Nt, 



(114) 
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4.4.2 The four point vertex function 

The expression of I3x can be obtained by considering the renormalization of the vertex 
function at first order in x and At this order, the four point vertex function is given 

using the diagrams of Fig. 8 by: 



r[fl/Ar (^1) ^2, ^3, ki) 



(^1 + k2)^ ih + 



^ Q d-2 r d'^k T^\k) 



5) 



The renormalized vertex function is defined through: 



^^R (ku k2, k3, k^) = rjf//^ (fci, /C2, h, h 



(116) 



where Z is the spinon field renormahzation. At this order we just need the field renor- 
malization Z of the 0{2N) model. Defining Z^t by: 



{xt)ii = Z^t xt 



(117) 



we deduce from eq. (|115|) : 



4(rf-l) f d<^k Tl^{k) 



Nd J {2'KY{k^ + w?y 



and explicitly: 



■'xt 



1 , A2 
1 - Tflnln — 

A m"^ 



1 8(c^-l) ^ A 
JSd m 



m d = 2 



in d> 2 . 



118) 



(119) 



Using these results, together with the expressions of the (3 functions for t given by 
eqs.( |113| , p,14]) , we deduce the /3 function for the variable x: 



xt f ^ t 



Nd 



t 



X 



in d = 2 



in d > 2 



(120) 



where in the last equation, (3^^'^^'^ is the (3 function for t for the 0{2N) model. 
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4.4.3 Comparison between weak coupling and large analysis 

In dimension two, one sees that tlie /? functions of eqs.( |113| , |12(J| ) agree, at first order in 



X 



and witli tlie weak coupling expressions of eg . (11061) . There is still a disagreement for 
the two loop term of the f3 function for x but this corresponds to subdominant corrections 
to the scaling behavior at the 0{2N) fixed point. Anyway, this term does not contribute 
to the mass gap at this order. Indeed, one can check that the expression of the correlation 
length obtained from eq.( p.l2| ): 

identifies with the correlation length (eg. (^5])) obtained in the weak coupling expansion. 
In dimension greater than two, the cross-over exponent (f)o{2N) is obtained from eg. (|120|) : 

0o(27V) = d-2 + ^^^Sd . (122) 

For example, in d = 3, we find (po{2N) = 1 + 32/37r^iV, in agreement with the result of 
Chubukov et 

alii El. 

As seen, 0o(2Af) is not given by dimensional analysis but displays 
a non trivial correction at order 1/A^ contrarily to what was found in the preceeding 
subsection. This is the result of the important fact that, in the 0{2N) regime, the gauge 
field mass does renormalize contrarily to what asserted by Campostrini and Rossi 
Finally, it is easy to show that eg. ( |122| ) identifies, at order 1/A^ and with the weak 
coupling result (eg.(|5^)): 

e 

(poi2N)=e + — + — . (123) 

Altogether these results show that the double expansion in 1/A^ and x agrees with the 
weak coupling analysis provided one is sufficiently closed to the 0{2N) fixed point. 



5 CONCLUSION 

In this work we have studied the SU{N) ® U{1)/SU{N - 1) O f/(l) NLa model which 
is relevant for the low energy physics of frustrated spin systems. We have shown that 
this model interpolates between two different fixed points which are 0{2N) and CP^^^ 
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symmetric. Our main goal, in this paper, was to inquire the rehabihty of perturbation 
theory in such a multi- fixed point situation about which httle is known at present. To 
this end, we have performed a rather complete analysis of the RG properties of the model 
using the weak coupling and expansions at next to leading order. 

The main characteristic of frustrated systems is the presence of an additional length 
scale S that determines the scahng regimes governed by the different fixed points of the 
theory. While the presence of this scale does not alter the weak couphng expansion, it has 
strong consequenses for the expansion. Indeed, the standard approach of the 
expansion fails to explore all the phase diagram: the decoupling procedure used to tackle 
with the U{1) current-current interaction is meaningless as one removes from the basin of 
attraction of the CP^~^ fixed point. We have shown that another strategy has to be used 
to explore the physics for large values of the gauge field mass M^: around the 0{2N) fixed 
point a double expansion in 1 /N and x is necessary to obtain the correct renormalization of 
the parameters of the massive CP^~^ model. Doing this, the results obtained from weak 
coupling and expansions agree in the very neighborhood of the two fixed points of the 
theory, giving us confidence on perturbation theory in the vicinity of these fixed points. 
However this procedure also addresses the question of universality for frustrated systems 
and for models where several antagonistic fixed points coexist. Indeed, the weak coupling 
expansion of NLcr models, as we know from the work of Fricdan, is renormalizable. This 
implies that a priori only two coupling constants x and T are sufficient to determine 
the low energy behavior when <^ p <^ A, i.e. in the whole range of variation of x, 
< X < 1. In this sense we expect the behavior of frustrated systems to be universal. 
This scheme appears to be more questionable in view of the large N expansion since one 
may suspect that in the double expansion in 1/N and x, higher orders corrections in 
X involve more and more irrelevant operators and coupling constants. If it is the case 
the existence of a universal field theory describing the whole phase diagram of frustrated 
systems should be ruled out and only a local idea of universality, bounded to the vicinity 
of the fixed points, would be meaningful. 

A pratical consequence of the previous discussion is that it could be very difficult to 
observe the 0{2N) scaling behavior when dealing with numerical simulations with systems 



39 



of size L <^ S. Indeed, for sizes of order L ~ S one can expect to observe rather than 
the 0{2N) scahng some effective exponents that result from the competition of 0{2N) 
and CP^~^ fluctuations. This might be one reason why, up to now, the 0(4) critical 
behavior has not been observed in the Monte Carlo simulations of frustrated Heisenberg 
spin systems. A better understanding of finite size scaling in the model is in progress. 
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Appendix : Gauge field renormalization. 

We give here the asymptotic UV behavior of the diagrams contributing to the renor- 
mahzation of the transversal part of the inverse gauge field propagator r±{p) in 2 and 
3 dimensions. Since there is no renormalization of the longitudinal part of the gauge 
field propagator we can forget, in the diagrams of Fig. 5 involving internal loop of the 
gauge field, the longitudinal part Si\\ and keep only their transversal part Si±. More- 
over to determine the renormalization of the propagator, we can restrict to search for 
the counter-term proportional to since it defines completely the renormalization of the 
propagator itself. 

• In 0? = 2, we have obtained: 



^(^—\n\n — + ^—\i— 



^ ^ Y247rm|j ^ rnj^ iSnrnj^^mj^ 247rm|j ^ ^ 



S4±^0 

^-^^^(24^^K^-^^-^ + ^^S) 

/ 1 Af ^ / A^ \ 1 / A^ 
Se±=p^- 2-ln— -— ^In 27rM2-2 + ln— +— ^ In 2^ - 2 + In — 
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The global contribution to ^±{p) at order 1/N thus writes: 

Si + 2S2 + 2Ss + S^i_ + 2S5± + 2Sei_ + iS7± = 

and therefore we have: 

r±i/iv(p) = r^p) - Y2^^ (i^in ^ + ~ ^""^'^ (^''^' - 2 + In 

• In d = 3, the cut-off dependence of the diagrams writes: 

^ ^ \67r^m|j 487r^m|j 37r^mfl 

g^^^2 ( A 1 ^^A^ , 1 

\ Qn^mji ASn^mji ^n^mR m|j 127r^mR rnjij 



'4± 



/ 4 . /A ,^2. 



and the complete contribution to rj_(p) at order is: 
Si + 2S2 + 2S3 + S^i_ + 2S5± + 2Sei_ + 4S7i_ = 



( + In ^ + -^{l&mR - 487rAM2) In + m^A 
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so that: 

/X ^ / X 1 / 32 , A 16A 256 /1 nAM^\ ^ /A , m 

' ASTrniRN \67r^ rriR TrniR tt^ \6 rriR J \lo ztt 

(129) 
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Figure Captions 



Fig.l: The infrared RG flow in ci = 2 in the {x,T) plane. 
Fig. 2: The infrared RG flow in o? = 2 + e in the (x, T) plane. 

Fig. 3: Feynman rules for the 1/N expansion in the CP^~^ regime. The solid line 
denotes spinon propagator while dashed and wavy lines represent respectively the A and 

Afj, propagators. 

Fig. 4: Diagrams contributing to the two point function at order 1/N in the CP^~^ 
regime. 

Fig. 5: Diagrams contributing to the gauge field propagator at order in the CP^~^ 
regime. 

Fig. 6: The four point vertex occuring in the double expansion in x and 1/N in the 
0{2N) regime. 

Fig. 7: Diagrams contributing to the two point function at first order in and x in 
the 0{2N) regime. 

Fig. 8: Diagrams contributing to the four point vertex function at first order in 1/N 
and X in the 0{2N) regime. 
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Fig. 1. 
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Fig. 2 




Fig. 3. 
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Fig. 4. 
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Fig. 8. 
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